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Abstract. We study possible effects of Lorentz-violation on the generation of 
cosmological perturbations at inflation by introducing a simple inflating five- 
dimensional brane-world setup with violation of four-dimensional Lorentz-invariance 
at an energy scale k. We consider massless scalar field, meant to mimic perturbations 
of inflaton and/or gravitational field, in this background. At three-momenta below 
k, there exists a zero mode localized on the brane, whose behaviour coincides with 
that in four-dimensional theory. On the contrary, at three-momenta above fc, the 
localized mode is absent and physics is entirely five-dimensional. As three-momenta 
get redshifted, more modes get localized on the brane, the phenomenon analogous to 
"mode generation" . We find that for k H, where H is the inflationary Hubble scale, 
the spectrum of perturbations coincides with that in four-dimensional theory. For 
k < H and time-dependent bulk parameters, the spectrum deviates, possibly strongly, 
from the flat spectrum even for pure de Sitter inflation. 
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1. Introduction and summary 

The standard inflationary mechanism of the generation of cosmological perturbations 
assumes Lorentz-invariance. This assumption is worth questioning [TJ EJ El HI El M El 
El 01 EDI EH [121 CHI US CHI CEl CZl EHl EHl EOl EH |22j especially in view of the fact that 
relevant perturbations initially have wavelengths much shorter than the Planck length 
(this observation is often referred to as "trans-Planckian problem"). Furthermore, there 
is no guarantee that Lorentz-invariance holds even towards the end of inflation. Thus, 
it is certainly of interest to understand the effects that Lorentz-violation may have on 
the cosmological perturbations. 

Several ways to address this issue have been proposed. One is to assume Lorentz- 
violating dispersion relation at short enough wavelengths [TJ EH EJ EJ EJ E]- If at 
asymptotically high momenta the frequencies exceed the inflationary Hubble parameter, 
the adiabatic vacuum is a natural choice of the initial state, and then the primordial 
spectrum is still flat for pure de Sitter inflation [TJEJEJIIIE], while the predictions for the 
tilt may be different, as compared to the standard mechanism, if the expansion is not 
exactly de Sitter. The opposite case of low frequencies at high momenta pjj [2] is hard 
to treat: perturbations are not initially in the adiabatic regime, so there is no preferred 
initial state (an exception here is the bouncing Universe model of reference jS]). 

Another approach is to make ad hoc assumptions on the initial state of 
perturbations [HI El UH1 EJJ • The problem with this approach is that any choice other 
than the adiabatic vacuum is hard to justify, whereas different choices lead to quite 
different spectra of perturbations [TUJ ITT] . 

One more possibility has to do with "mode generation": there are models where 
the number of modes increases in time, as the Universe expands [T2* ] ll3 | IT H ITo 1 ll6j: this 
happens, e.g., if the physical momentum is effectively bounded from above. This option 
is in a sense similar to the previous one: the result for the spectrum of perturbations 
strongly depends on the choice of a state in which a newborn mode emerges. If this is 
the adiabatic vacuum, the outcome is similar to the standard scenario, at least for pure 
de Sitter inflation [3 d EH] • 

The emerging picture is that in the least exotic scenarios, the spectrum of 
perturbations is flat for pure de Sitter inflation, while predictions for the tilt in non- 
de Sitter cases are generally different from those of the standard mechanism. If 
the energy scale k associated with Lorentz-violation is much higher than the Hubble 
parameter H at inflation, the deviations of the spectrum from the predictions of the 
standard mechanism are suppressed by powers of H/k, as argued on general grounds, 
e.g., in references jH EH EE]. O n the other hand, dramatic assumptions, like, e.g., 
the Corley-Jacobson dispersion relation [2B] with frequencies tending to zero at high 
momenta, lead to dramatic consequences: the spectrum may be quite different from 
the flat spectrum even for pure de Sitter inflation [TJ E]. Also, for k < H the 
properties of the primordial perturbations may deviate substantially from the standard 
predictions even in some less exotic models, like, e.g., the model of reference ^H] based 
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on noncommutative field theory. 

Violation of four-dimensional Lorentz-invariance is fairly natural in brane-world 
models [2H E51 ESI EH EH] • It may occur in the bulk, especially if there is bulk matter 
(a particularly well-understood solution of this sort is a brane moving in Schwarzschild- 
anti-de Sitter bulk [201120]) or Lorentz- violating bulk condensate pTT] . Since the wave 
functions of brane matter may extend to the bulk, the dispersion relations on the brane 
are then naturally Lorentz- violating |2Zj. Furthermore, some brane-world geometries 
have the property of "mode generation" in the following sense: for low three-momenta, 
there exist modes localized (or quasi-localized) on the brane, while for high three- 
momenta, localized and even quasi- localized modes are absent |2HJ - Thus, as the 
Universe expands and three-momenta get redshifted, more and more brane modes 
emerge; this may be viewed as the brane-world version of mode generation. Yet the 
initial state is well defined, provided that modes with high three-momenta are adiabatic 
in the bulk. 

In this paper we consider an inflating version of the setup of reference |2E] , and study 
massless scalar field, meant to mimic perturbations of the inflaton and/or gravitational 
field, in this background. This is nothing more than a toy model, as we specify neither 
the mechanism giving rise to Lorentz-violation, nor the mechanism driving inflation. In 
fact, our background may be totally unrealistic, as it may require matter with unphysical 
properties either in the bulk or on the brane, or both. Also, since we do not have a 
mechanism generating our background, we cannot quantitatively address the issue of 
back reaction (201 HI EH E2] of the perturbations on background geometry (however, we 
will argue on qualitative grounds that this back reaction is small). On the other hand, 
the advantage of our approach is that the properties of the massless scalar field are well 
understood in spite of mode generation on the brane and other unusual features. 

To a certain extent we confirm the picture emerging from the previous studies. In 
the simplest version of our setup, with time-independent bulk parameters, the spectrum 
of perturbations generated due to inflation is flat for pure de Sitter expansion, while its 
amplitude and, in the case of non-de Sitter expansion, the tilt are generally different 
from those predicted by the standard mechanism. Furthermore, for k ^> H, where k and 
H are still the scales associated with Lorentz-violation and expansion rate at inflation, 
respectively, even the latter deviations are small. 

Things are more interesting if one makes a relatively mild assumption that the bulk 
parameters depend on time, so that k is smaller than H towards the end of inflation. 
As an example, we consider the case in which k slowly decreases in time, from k 3> H 
at the early stage of inflation to k < H towards its end. Then the primordial spectrum 
is flat at low momenta, gets strongly tilted at intermediate momenta, and flattens out 
again at high momenta, as sketched in figure ^ The position of the tilted part, as 
well as its width, depend on the parameters of the model; this part may well extend 
over several (but not necessarily many) decimal places, and the amplitudes at low and 
high momenta may differ also by several orders of magnitude. To avoid confusion, the 
primordial spectrum is cut off in the ultraviolet in the usual way, the cutoff being at 
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physical momenta of order H at the end of inflation. 

The opposite situation, shown in figure 121 occurs if k increases in time from k < H 
to k > H. 

A lesson we learn from our toy model is twofold. On the one hand, we confirm that 
it is quite unlikely that in a well defined setup with Lorentz-violation well above the 
inflationary Hubble scale (k 3> H in our notations), the standard predictions may be 
strongly modified. In our model, physics at three-momenta above k is entirely multi- 
dimensional, and yet fork^>H this has basically no effect on primordial perturbations. 




Figure 1. Power of perturbations, as function of momentum p, generated in a model 
with k decreasing in time. p x is a free parameter and e <C 1 is another free parameter 
of the model. The exponent a is model dependent and varies in the range < a < 3/2. 
See section T5. II for explanation and qualifications. 




Figure 
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On the other hand, if Lorentz-violation is characterized by low energy scale k < H, 
the spectrum of perturbations strongly depends on details of the cosmological evolution 
at inflation, and may well be grossly different from the standard approximately flat 
spectrum. 

The spectrum sketched in figure Q is interesting. If our scalar field is interpreted 
as mimicking tensor perturbations, this spectrum corresponds to greater (maybe much 
greater) primordial amplitudes of gravity waves at shorter wavelengths^. If so, a possible 
non-observation of the effects of tensor perturbations on CMB would not preclude the 
detection of primordial gravity waves at shorter wavelengths by techniques ranging from 
timing of pulsars to ground-based interferometers. Moreover, since the latter techniques 
cover a wide range of wavelengths, one might be able, with luck, to uncover an interesting 
pattern of relic stochastic gravity waves in our Universe. 

Likewise, the spectrum shown in figure El would correspond to relatively high power 
at long wavelengths, possibly measurable through CMB, with deviations from flatness 
potentially detectable by experiments searching for gravity waves at shorter wavelengths. 

We conclude this section by mentioning that even more complicated spectra are 
obtained if one allows for time-dependence not only of k, but also of other bulk 
parameters. We will comment on this point in section 



2. Setup 

We consider (4+l)-dimensional model and denote the coordinates as (rj, x\ y), i = 1, 2, 3. 
We take the five-dimensional metric in the following form (signature +, — , — , — , — ), 

ds 2 = a 2 (7]) [a 2 (y)dr] 2 - f3 2 (y)dx 2 } - a 2 (y)dy 2 (1) 

where we made a convenient choice of the coordinate y. There is a single brane at 

y = Vb = o 

As usual, the warp factors a(y) and (3{y) are continuous across the brane, while their 
derivatives d y a and d y (3 are not. For equal warp factors, a(y) = /3(y), and static 
geometry a(r]) = const, equation (0) has a general form of warped metric invariant 
under four- dimensional Lorentz group. For a(y) ^ P(y), and still airf) = const, the 
coordinates can be chosen in such a way that the induced metric on the brane is still 
four-dimensionally Lorentz-invariant, while the full five-dimensional metric is not. This 
is the case studied, e.g., in references (23 EB]. in what follows we choose coordinates 
so that the brane is Lorentz-invariant for a = const. To specify our setup further, we 
assume that both a(y) and (3(y) are Z 2 -symmetric across the brane and decay away 
from the brane, 

<x(y), P(y) , asy^oo 
and that their ratio tends to a small constant, 

W) ^e, as^-oo (2) 
| As usual, the gravity waves decay as a -1 after re-entering the horizon. 
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The case e = and a(r]) = const was discussed in reference (2HJ- 

We are interested in inflating setup. For de Sitter inflation, which is the main focus 
in this paper, one has 

«*!) = ~ (3) 

If inflation is not exactly de Sitter, the Hubble parameter H slowly varies in time. Later 
on we will also consider generalizations of the setup (|T} in which a and/or (3 slowly vary 
in time as well. 

Throughout this paper we concentrate on the properties of five-dimensional massless 
minimal scalar field <fi in the background (0). Its action is 



S<t> = \j d 5 X y/gg AB d A cl>d B < 



Upon introducing another field \ via 
1 



r a(3 3 / 2 
the action becomes 



X (4) 



a 



(xT + ( % + % ) x 2 



S x = -J dnd'xdy ^ + -X' - ^(«dt) J - r , ■ V2 ■ 4 p 

where dot and prime denote d v and d y , respectively. The field equation in terms of 
three-dimensional Fourier harmonics is 

-X+-X- U(y)p 2 X ~ a 2 ( V ) [- X " + V(y) X ] = (5) 

(X 

where p is conformal three-momentum, and the "potentials" are 

u(v) = m 

and 

3(3" 3 P' 2 
Viy) = 2j + AW 

In what follows we assume that V(y) tends to a constant as \y\ — > oo, which, for future 
convenience, we denote as 

v(lvl - °°) = 

is a monotonically decreasing function whose asymptotics, according to 
equation (J2J), is 

oo) = e 2 

One more assumption is that the width t/ of the "potential" {/(?/), i.e., the size of the 
region where U(y) substantially exceeds e 2 , is sufficiently small; roughly speaking, 

yo < k- 1 (6) 
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The reason for this assumption will become clear later§. Finally, we will assume that 
k and e are the only relevant parameters entering U(y) and V(y), the width y being 
roughly of order k^ 1 or maybe substantially smaller. 

We have to specify the boundary condition on the brane. We impose Z 2 symmetry 
across the brane, consider the theory on half-space y > 0, and, guided by the 
properties of gravitational perturbations in the Randall-Sundrum model (see, e.g., 
and references therein), impose the condition 

<j>'(y = 0) = 



X'-~X) =0 (7) 



This is in accord with the discontinuity of /3'(y) at y — 0, which results in the 5-function 
behaviour of V(y), 

V(y) « 3^5(y) at ^0 (8) 

Hereafter we assume that /?' < for all y, so that this contribution to the potential is 
attractive. 

Although our analysis will not depend on the precise form of the warp factors a(y) 
and P(y), in section |H1 we will give an explicit example by choosing 

m = e~ ky (9) 

Then 

9 

V(y) = -k 2 = const 

Equation (|SJ) is simplified further for a particular choice 

U(y) = ^S(y)+e 2 (10) 

(the reason for the coefficient 2/ (3k) here will become clear later). With this choice, the 
field equation and boundary condition read 



2 9 2 / \ 

X + -X - e V X ~ a iv) 

a 



-X" + \k 2 x 



(11) 



, 3, 1 p 2 \ 
2 3ka 2 { V ) J y=Q 

Note that even though variables separate in the field equation (|llj). the dynamics is 
non-trivial because of the time-dependence of the boundary condition ([12)1 . 

§ We are grateful to W. Unruh who pointed out the relevance of this assumption, which simplifies the 
analysis considerably. We expect, however, that our main results remain valid at yo > and even 
yo S> We do not attempt to perform the anaysis in the latter cases, as our attitude in this paper 
is "proof by example" . 
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In static background, a{rj) = 1, and with \ — e~ tuJV Xu> equation (J5J) and boundary 
condition (JJj) take the form of the static Schrodinger equation 

- x!L + W p {y)xu = ^X. (13) 

where x{v) m ust be symmetric across the brane, and 

W p (y) = V{y) + U(y)p 2 

(the 5-function term (jHJ) is implicit in V(y)). This effective potential has volcano shape, 
as shown in figure El At p — 0, continuum of bulk modes starts at to = and there is 
a bound state with u = whose wave function is 

Xo(y) = NP*(y) (14) 

where N is a normalization factor, ensuring 

oo 

j \xo{y)\ 2 dy = l 
o 

This solution corresponds to (f>o(y) = const. We assume that this is the only bound 
state at p = 0. 




Figure 3. Volcano-shaped effective potential W p . Its height and large- y asymptotics 
increase with p, while the dip at y — is independent of p. 

By continuity, at p k the bound state is still present. Depending on the shapes 
of U(y) and V(y), there may or may not exist quasi-localized modes at p ~ k. For 
sufficiently small width of the "potential" U(y) (see ©), at p 3> k even quasi-localized 
modes cease to exist, while continuum of bulk modes starts at u 2 = e 2 p 2 + \k 2 . These 
properties are illustrated in Appendix. 
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As an example, the localized solution to equations and (J 12)) has the form 
Xo = iVexp ^ [^--j \ 
and exists only for 

3 , 

At p -C k the dispersion relation is u 2 = p 2 [l + 0(e 2 )} which explains the coefficient in 
equation (fTfjjh 

4. Inflation with time-independent parameters 

4-1. Preliminaries 

At large negative times, the third term in equation (jSJ) dominates, so all modes are 
shifted in the bulk towards large y, where U(y) = e 2 . For given p, all frequencies tend 
to u = ep as rj — > — oo. Thus, at large negative times the evolution proceeds in the 
bulk and is adiabatic; the field \ has well defined decomposition into positive- and 
negative-frequency components, 

d\ 

where oj\ p {rj) and X\,p(y> v) are eigenvalues and properly normalized eigenstates of the 
operator 



X 



Xx 



j,exp ^-i J uJx, P dr]j A Xp + h.c. 



(15) 



U(y)p 2 + a 2 ( V )[-d 2 + V(y)} 



(16) 



and A^ p are creation and annihilation operators. The asymptotics of xx,p at large y is 

Xx, P -> cos(Ay + (p x , P ) 
where the phase (p\ tP depends on r/. Thus, 



9 



< P (v) = a 2 (v)[X 2 +^-k 2 )+e 2 p 2 



(17) 



As the physical momenta 



P 



a(r}) 



get redshifted, the properties of modes change. At late times the third term in 
equation © is negligible for superhorizon modes (P <S H), and besides the bulk modes, 
there is a localized mode 

xioc(y,v) = N{3 3 Hy)Mv) (18) 

where ipoiv) obeys the usual equation for superhorizon modes in four dimensions, 

a 

-ipo + -ip = 
a 
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and has a growing component 

Vto = 0o • 0(77) (19) 

where 0o is independent of time. This is the standard superhorizon mode. It is localized 
on the brane, the wave function in the transverse dimension being given by equation ()14j) . 
In terms of the original field 0, we have for the localized part 

01oc = 00 

Our goal is to relate the constant O i n equation (fT9"j) to the creation and annihilation 
operators entering equation (|T5|) . This analysis is different for k S> H and k < H . 



4.2. A; > H 

It is straightforward to see that in the case k 3> H, the evolution of the scalar field is 
adiabatic when the physical momentum P(rj) becomes of order k, since at that time the 
physical frequencies are high, 

and slowly varying, 

< H 



an 2 ~ k 

At the time when P ~ k, a bound state with f2 ~ ^> H emerges in the spectrum 
of the operator (|16|h and then its frequency o>(?7) slowly decreases from u ~ ak at the 
time when P ~ /c to = p at later times when P k but still P ^> H. In this 
time interval the evolution is still adiabatic. Since there is a gap of order k between 
positive- and negative-frequency parts of the spectrum, there is no mixing between the 
positive- and negative-frequency modes at the time when the bound state appears, so 
the bound state is initially in the adiabatic vacuum, up to corrections suppressed by 
H/k. Therefore, its further evolution proceeds in the same way as in the standard four- 
dimensional scenario; at the horizon crossing, P ~ H, the shape of the bound state is 
already given by equation (|TH|) modulo if //c-corrections. It is important here that by 
the time of the horizon crossing, there is no cross talk between the bound state and 
bulk modes, since they are also separated by the gap of order k (in terms of physical 
frequencies) . 

We conclude that for k S> H, the spectrum and amplitude of perturbations coincide 
with the standard predictions up to corrections suppressed by H/k. 



4.3. k<H 

In this case the physical momentum crosses H at the same time as it crosses k, or earlier. 
Adiabaticity gets violated at the same time or before the bound state appears. This 
case is hard to treat in detail (see, however, sections El and EJ) , but for pure de Sitter 
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inflation, equation (j3J), there is a simple scaling argument showing that the spectrum in 
the localized mode is flat. Let us introduce a new time variable 

C=PV 

Equation (jHJ) written in terms of ( does not contain p explicitly. The initial condition is 

x K _L e -^ [ d\( X x(y, C)A- Xp + h.c.) 

where the factor (qo) -1 / 2 comes from (uj\, p )~ l l 2 in equation ()15j) . and we recalled that 
^Ajp — e P as V — > —oo. Thus, up to overall factor p" 1 ^ 2 , the initial condition does not 
contain p explicitly either. The boundary condition (0) is independent of p. So, the 
solution to equation (JHJ) has the general form 

* = > c) 

where for given £, the operator F is a linear combination of Af p with p-independent 
coefficients. Therefore, towards the end of inflation, the localized part of the solution, 
for sub horizon three- momenta, behaves as follows (see equation (JTHj) ). 

xioc = *Po(v)xo(y) 

where Xo(y) is given by equation (jH|) and 

Mv) = ^io p 
' id. 



^3/2 ^ 

where O p is a time-independent linear combination of Af p with p-independent 
coefficients. The latter property implies that O p describes a Gaussian field obeying 

(OpOp'> = ^ 3 (p-p') 

where C is independent of p. Recalling the relation (fl|). we obtain that the localized 
modes 

0o(p) = ^o(p)a _1 (^) 

towards the end of inflation make up a Gaussian field with flat spectrum, 

const 
p 3 

It is clear that this relation holds only for momenta which are superhorizon by the 
end of inflation, higher momenta being suppressed in the usual way. Therefore, if the 
overall amplitude of perturbations is much smaller than 1, their back reaction on the 
background geometry is likely to be small. The latter remark applies also to more 
interesting cases studied in section 03 

It is also clear that for k < H, the constant in equation (J2U|) depends on k and 
H in a non-trivial way (see also sections El and EJ). The latter observation implies, in 
particular, that for non-de Sitter inflation, the expression for the tilt does not generaly 
coincide with that in the standard scenario. 



>o(p)0o(p')> = ^5 3 (P - p') (20) 
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5. Inflation with slowly varying k 

5.1. Decaying k 

Let us now consider inflation with time- dependent bulk parameters. As an example, let 
us assume that k slowly decreases, while, for definiteness, H stays constant at inflation. 
Let k cross H at some rj = n x , as shown in figure EJ We assume in what follows that 
k is not grossly different from H after t] x , and will not count factors of H/k. On the 
other hand, we consider e as a small parameter, and will be interested in the dependence 
of perturbations on this parameter. The analysis of modes with low and high three- 
momenta gives different results. 




Figure 4. Decaying k. 



(i) Low three-momenta. These are modes that cross out the de Sitter horizon 
before the time rj x , i.e., 

P(r,x) = ~rr <: H (21) 

The analysis of these modes goes through as in section 14.21 the spectrum is flat and 
the amplitude of perturbations in the localized mode is given by the standard four- 
dimensional formula 

MP) ~ ^ (22) 

(ii) High three-momenta. These modes obey 

eP(77x) > H (23) 

We note in passing that there is an intermediate range of three-momenta, such that 
H <^ P{j] x ) <^ He^ 1 . This range is actually wide, since we assume e < 1. 
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High-momentum modes (|23|) are adiabatic at rj = i] x . We will see that the 
interesting range of transverse momenta is A <C k. Then the interesting modes get 
out of the adiabatic regime at the time i] e when 

eP(Ve) ~ H (24) 

(we call this moment "e-horizon crossing" ; note that this moment occurs after k crosses 
H, i.e., 7] t > r]y). At this time P(rj e ) 3> k, so for A < k the modes are still far away from 
the brane, and the frequencies are all of order ep. Just before the e-horizon crossing, 
the frequencies (fTTjl are 




\2 I 9p 

u x , p = ep\l 1 + R i (25) 



where we made use of equation (}24j) to write 
a (Ve) = Jj 

In what follows we will often omit the second factor in equation ()25|h since it is of 
order 1. 

The frequency (|23j) is of order a(rj t )H, which means that adiabaticity indeed just 
becomes violated. According to equation (|15jl . the amplitude of x a t this moment for 
all interesting A is 

X(Ve) oc 
oc 



1 



(ep) 1 / 2 

During some time after e-horizon crossing, the modes are still away from the brane, so 
equation (jJJJ) may be approximated as (modulo mixing|| between modes with different 
A) 

.. 2 1 A 2 + jfc 2 
~Xx + ~^X\ o TH — ^ A = 

For A <C k (and k < H, recall that r] e > r/ x ), this equation has both decaying and growing 
solutions. The growing solution, which is of primary interest, behaves as follows, 

Xx oc v - a ^ +l ' 2 (26) 

with 



a(X) 



A 2 + Ik 2 



H 2 

This behaviour terminates at the time rjk when 

P(Vk) = -j^r ~ k (27) 

|| This mixing can be shown not to change the estimates that follow, provided that fc is not grossly 
different from H, which we assume throughout. 
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Around this time the wave functions change considerably; in particular, the localized 
mode emerges. 

The amplitude of x a t this moment is 

, X«(A)-l/2 

Xx(Vk) ~ x(Ve) — 
\VkJ 

H \ q(A)-1/2 



(ep)V2 \ek^ 

With e C 1, the integral over the transverse momenta A is saturated at A <C k, which 
justifies the assumption made earlier. Since all relevant parameters (P, k and H) are 
roughly of the same order of magnitude at this time, the estimate of the amplitude of 
the localized mode at rj ~ rjk is, up to logarithms, 

where Xo(y) is still given by equation (JHj), 



3 / k 2 

a = a(A = o) = -yi--p 

and f(H/k) is roughly of order 1. Let us note that for very small three-momenta, the 
value of k entering here is the asymptotic value at rj ^> i] x . 

Later on, at r/ > r] k , the localized mode grows as a(rj) (since k < H, this mode is 
superhorizon at rjk and later). Hence, towards the end of inflation the amplitude of the 
localized mode is (we use the notation (fT£|) ) 

, , 11 a( V ) (H 
Mri) ~ ^pT^) S U 
Recalling equation (|27j) . we obtain in terms of the original field <fi that the localized 
mode of perturbations stays frozen, and its amplitude is 

Mp) = ^f(j) (28) 



where / is again roughly of order 1. We obtain that at high momenta obeying the 
relation ([23)1 . the spectrum is again flat, while the amplitude is enhanced by the factor 

e~ a (29) 

as compared to the modes in the low momentum region, equation (|22|). Note that 
< a < 3/2 with a -> 3/2 as fc/i? -> 0. 

The overall behaviour of perturbations is shown in figure ^ Let us make several 
remarks concerning this result. First, the enhancement of the amplitude, as given 
by ([29)1 . is a free parameter, and can be of several orders of magnitude. The reason for 
this enhancement is clear: modes with low three-momentum become "tachyonic" when 
they are still living in the bulk. For e = and constant k < H, all modes are tachyonic 
"from the beginning" , the situation reminiscent of that in four-dimensional theory with 
the Corley-Jacobson dispersion relation (HEj- With our setup, we effectively cut off 
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A 




Figure 5. Increasing k. 



the divergence associated with the tachyonic instability. Second, the beginning of the 
transition region p x is determined by the time rj x when k crosses H. Clearly, the 
physically interesting situation occurs when this time is not very distant from the time 
inflation ends. Third, the momenta in the end and in the beginning of the transition 
region differ at least by a factor e _1 (this is the case shown in figure^), but their ratio 
may be much larger, since it is not necessarily determined by the parameter e only, as 
one might guess by comparing the relations (j21j) and (J2"3j) . The value of the parameter 
k entering the exponent a is to be evaluated somewhere in the interval (r] € (p), f]k{p)), 
so for very slowly varying k(r)) this value in fact depends on three- momentum p and 
settles down only as k(r]k(p)) reaches its asymptotic value. Fourth, the parameter e 
may itself depend on time, which would lead to even more complicated spectrum of 
perturbations. Finally, although we have considered an example with k crossing H, an 
interesting spectrum is obtained also if k < H during entire inflationary epoch, provided 
that k and/or e depend on time. The latter two points are clear from equation (J2*5|l . 
which is valid for k < H, and the above observation that e and k entering equation (f2*%|) 
depend on three momenta, as they are to be evaluated at some r\ belonging to the 
interval (rjeip), Vkip))- 



5.2. Increasing k 

If k increases in time as shown in figure the properties of soft and hard modes (in 
three-dimensional sense) are reversed. Soft modes in this case are those for which the 
bound state gets formed before k crosses H, i.e. 

It is these modes that become tachyonic when still living in the bulk. They increase 
as given by equation (|26|) from the time i] t to r]k', these times are still determined by 
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equations (J23J) and ()27|). with k now being the value of k{rf) in the past, where k < H. 
At the time % the bound state emerges, and then the amplitude of the bound state 
mode increases as aiji) -1 . Thus, equation (|2*%|) applies now to the soft modes. 

Hard modes are those which never become tachyonic when staying in the bulk, 
which is the case for 

eP(7? x ) 

They start growing after forming the bound state, so their amplitudes are given by the 
four-dimensional expression fj22(l . The overall spectrum is thus the reverse of figure P 
it is shown in figure El The remarks made in the end of section 15.11 apply to the case of 
increasing k as well. 



6. Explicit example 

In this section we consider an explicit setup outlined in the end of section 121 The 
field x obeys equation (JTTJ) with the boundary condition (JT5J). We will be interested 
in the case k ~ H and, as in section study the dependence of the amplitudes on e, 
assuming e < 1. As mentioned in section |21 in spite of the fact that variables separate in 
equation (JTTJ), the dynamics is non-trivial due to the time-dependence of the boundary 
condition ()12j) . We simplify the analysis further by making use of the instantaneous 
approximation to the boundary condition ()12j) . Namely, at rj < where 7] k is given by 
equation (|2"Tj) . i.e. 

Vk = ~ 1 (30) 

Hp p 

the third term in equation (|12|) dominates and we omit the first two terms. So, at these 
times the boundary condition reads 

x\ y =o = at r) < r] k (31) 

On the other hand, at r\ > T]k we neglect the third term in equation (fT2"|) and the 
boundary condition becomes 

x' + ^kx) =0 at r 1 >n k (32) 

1 J y=0 

In this approximation, the problem reduces to that of solving the field equation fTTf 
with the boundary conditions (f3*Tj) . (|3~2|) and matching the solutions at r] = r^. 

At r] < r] k the solution is (we slightly change the notation as compared to 
equation (fT5~|) ) 



oo 

x = / -7= [x\*{y> ^ A x, P + hx -] 

J V Z7T 



where x\,p are 

X\j>(v, X) = -iv^e^v^r/tf^-^ep) sin(Ay) (33) 
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and 

^fii'-ipj-iP < 34) 

As rj — ► — oo, these modes approach properly normalized negative- frequency plane waves 




Xx, P (y,v) = \^~^ + ^sm(Xy) 
V to 



where uj = ep and <p is a constant phase. For k < H and small enough A the parameter 
a is real. So, these modes indeed grow like |^| 1//2_a in the interval r/ e < 77 < where 77,= 
is determined by equation (|24jl . For k > H, the parameter a is imaginary, so the modes 
oscillate in 77 for all 77 < rjk. 

At late times, rj > rjk, besides the bulk modes, there are two localized modes obeying 
equations (fTTj) and (|32~j) . 

xtodviv) = Vsk ■ V^vi-v)^?'^ 

These are the standard four- dimensional modes; recall that k ~ H, so at 77 > 77^ these 
modes are superhorizon. At 77 > 77^ the field x is a superposition 

X(y, V) = C^XiociVi V) + C i+) xt oc (y, v) + bulk modes (35) 

The coefficients here are obtained by matching this solution to the solution (f3*3|) at 
V — Vk- We are interested in the coefficient C^~\ which is the amplitude of the growing 
localized mode. To this end we notice that the bulk modes entering equation ()35|) are 
orthogonal to Xio C with the standard scalar product 

= / dy(xip* - ip*x) 







On the other hand, 

(Xtoc: Xloc) = 3 

Thus 



where x is the solution ()33|). and the scalar product is to be evaluated at 77 = rjk. One 
obtains for arbitrary 77 



Thus 



00 

C(-) = / ^ e ^f -# {1) fe) - e<9 # (1) fe^ ^ A A~ +h c f36) 





where we took into account equation (J3U|). 
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Now we are ready to calculate the spectrum of perturbations in the localized mode, 
defined as 

2tt 2 
p 3 

We make use of the relation (@J) of the original field to the field x an d the standard 
commutation relations 



(p)<k,(p')> = — W(p-p') 



With the expression (J3l)j) for the coefficient we find 



4vr 2 7 3(A 2 + |A; 2 ) 4 



(37) 



For time-independent /c, the right hand side of equation (|3*7)l is independent of p, the 
spectrum is flat, while the amplitude depends on k, H and e in a non-trivial way. This 
is in accord with the general analysis given in section 0] 

On the other hand, if k slowly varies in time, the spectrum is no longer flat. The 
important time interval is T] e < rj < rjk, in which bulk modes grow provided that k < H, 
and do not grow for k > H . 

Let us consider the case of decreasing k, figure EJ For low three-momenta one 
has (see equation (J2T|) ) r] x ^> r] k and a(k) defined in (|3*l)l is pure imaginary. So, the 
integrand in equation (pTFJl is independent of e at small e. The spectrum is flat, and the 
power is 

H 2 

V (P) = 4^2 • const ' 

where the constant is of order 1. In the opposite case of high three- moment a, the times 
are ordered as r/ x < i] t < r] k . So, k(r] k ) < H and a(k, A) is real at small enough A. In 
this case 

2 1 ° W 4,r V«/ V 2/ T(3/2) 2 

Substituting this expression into equation (p?7j) one finds that the integral is saturated at 
small transverse momenta, A 2 ~ H 2 /\ne. Therefore, we again obtain the flat spectrum 
but with the amplitude enhanced by the factor (j^Hj) (up to logarithm), in agreement 
with the discussion in section HTT1 

If k increases in time as shown in figure a(k) is imaginary for high momenta 
(r] x > rjk) and real for low momenta (r] x < r) e ). The perturbations get enhanced by the 
factor e~ a at low momenta, and do not get enhanced at high momenta, so we obtain 
the spectrum shown in figure El 
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Appendix 

In this Appendix we illustrate the properties of eigenmodes of equation (JT3j) by 
considering a step-function "potential" 

U(y) = a 2 9(y -\y\) + e 2 9(\y\-y ) 

where a ~ 1 and e < 1. We take the warp factor f3(y) in the form ©, so that another 
"potential" is constant, 

v(y) = h 2 

We consider equation (|T3|) on half-line y > and impose the boundary condition (jJJ), 
i.e., 

X' + Ux) = (38) 

We are interested in localized or quasilocalized modes whose general expression at y < yo 
is 

= Ae~»y + Be»y , y < y (39) 

Then 

oo 2 = -v 2 + -k 2 + a 2 p 2 (40) 
while the boundary condition (}38|) reads 



3 k-u] A+ i-k + v ) B = (41) 



2 J \2 
Let us begin with localized modes. They decay at y > y , 

X. = e~ x( ~y-y^ , y > y 
where we set the overall factor equal to one. From equations (|T3j) and (jjQj) we obtain 

z/ 2 = A 2 + a 2 p 2 (42) 
where a 2 = a 2 — e 2 . By matching Xlo and at y = yo we find 

A = ^±V° 

b = ^ e -yo 

2v 

and the relation (|41|) gives 

— A u + 
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This equation, supplemented with equation (J42)) . determines the eigenvalue v for given 
three-momentum p. 

At p 2 = equation (fTH|) has a zero mode xo = exp(— \ky) with A = v = |fc and 
u) = 0. For p 2 <C A; 2 there exists a solution to equations (j42j) and (J4*3|) describing a 
localized state close to the zero mode, 

v = 3 -k + ^-e-^.p 2 + 0(p 4 ) 



2 - (a 2 - a 2 e~ 3kyo ) -p 2 + 0(p 4 



Up to rescaling of coordinates on the brane, this mode has relativistic dispersion relation. 

At p 2 3> k 2 no real solutions to equations (|4^|) . (f4*H|) exist: equation (j4^|l gives 
v > p, so the right hand side of equation (f4"Hj) is greater than 1 in that case. Thus, there 
are no localized modes at high three-momenta p. 

Let us now turn to quasilocalized modes. They still have the form at y < y$, 
while at y > they oscillate 

Xu = oe ii(y-y°) + * e -ig(y-yo) 

so that 

2 ~2 2 2 

g =a j) — v 

The normalization condition appropriate for the modes in continuum is 

\b? = ~ (44) 
By matching an d xL a ^ V = Vo we obtain 

1 1 u 

b = - (Ae~ uy ° + Be vyo ) + — (Ae"^ - Be uyo ) 

2 ' 2q v 7 

so the normalization condition (J4*4"j) and equation (J4*T]) determine the coefficient A, 



i v—\k \ 2 v 2 ( v— \k \ 2 
L4 2 = 2 e~ vyo + l-e vyo + — e~ uyo i-e^ 45 

The defining property of a quasilocalized mode is that it is large at the brane position 
(when normalized according to (|44|) ) . 

\A\ > 1 (46) 

At kyo ^> 1 this relation indeed holds for 

v = -k + 0{e~ 3kyo ) 

In that case the quasilocalized mode exists for arbitrarily high three-momenta p; 
furthermore, at p ^> k and 

v = -k- 3ke~ 3kyo + O —e~ 3kyo 
2 \p 



one has 



Ml ^ P e 3ky /2 

k 
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so the localization gets stronger at high momenta p. 

At kyo < 1 the only way the relation ()46|) may be satisfied is that p is large, and 
the first term in the square bracket in ()45j) vanishes, which gives 



c 



o Ik - v 



The latter equation, however, has positive solutions only for 

yo > h 

Thus, a quasilocalized mode exists at high three-momenta p provided that ky is large 
enough, otherwise there are no no quasilocalized modes at p ^> k. It is the latter case 
that is considered in the main text. 
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